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It is shown how suitably scaled, order-m moments, of the Elsasser vorticity helds in three- 
dimensional magnetohydrodynamics (MHD) can be used to identify three possible regimes for so¬ 
lutions of the MHD equations with magnetic Prandtl number Pm = 1- These vorticity fields are 
defined by = curl = to Tj, where are Elsasser variables, and where oj and j are, respec¬ 
tively, the fluid vorticity and current density. This study follows recent developments in the study of 
three-dimensional Navier-Stokes fluid turbulence [Gibbon et al. Nonlinearity 27, 2605 (2014)]. Our 
mathematical results are then compared with those from a variety of direct numerical simulations, 
which demonstrate that all solutions that have been investigated remain in only one of these regimes 
which has depleted nonlinearity. The exponents that characterize the inertial range power-law 
dependencies of the z’^ energy spectra, are then examined, and bounds are obtained. Com¬ 

ments are also made on : (a) the generalization of our results to the case Pm A 1 fl’Hd (b) the relation 
between Dl^ and the order-m moments of gradients of magnetohydrodynamic helds, which are used 
to characterize intermittency in turbulent hows. 

PACS numbers: 47.27.Ak, 52.30.Cv, 47.27.ek, 02.30.Jr 


I. INTRODUCTION 


which scales as 


Intermittency is widespread in nature: its char¬ 
acterization is a central problem in turbulence [ll- 
1^, nonequilibrium statistical mechanics, the pro¬ 
duction and storage of wind and solar energy, and 
the behaviors of market crashes, and of several crit¬ 
ical phenomena [13, mi- It has also been studied 
extensively in fluid turbulence P, PP and in mag¬ 


netohydrodynamic (MHD) turbulence |12h17| . of¬ 
ten by using order-p structure functions of fields 
such as the velocity and, in MHD, the magnetic 
field. An example is the (longitudinal) velocity (it) 
structure function, 


Sp{r) = ([(5M(r)]P), 

Su{r) = [u(x -hr) — u{x)] f , (1) 


* Postprint version of the manuscript published in Phys. 
Rev. E 93, 043104 (2016). 


Spir) 


( 2 ) 


for r]d r L, where iqy is the dissipation length 
scale below which viscous dissipation is significant, 
L is the large length scale at which energy is in¬ 
jected into the fluid, and the multiscaling expo¬ 
nents Cp, which are nonlinear, monotone increasing 
functions of p, characterize the intermittency [l|. 
Simple scaling is obtained if (^p depends linearly on 
p, as in the phenomenological approach (K41) of 
Kolmogorov [18|. 

To determine Cp is a challenging task H, [3 0 , 
which is especially difficult for time-dependent 
structure functions |l9|, l20| or MHD turbu¬ 


lence Therefore, we explore other signa¬ 

tures of intermittency. For three-dimensional (3D) 
fluid turbulence Refs. 21 have introduced a 
new way of analyzing direct numerical simulations 
(DNSs) to obtain fresh insights into suitably scaled 
(see below), order-m moments Dm of the vorticity 
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u> = V X u. These studies show the following: (a) 
on theoretical grounds, three regimes, I, II, and III, 
are possible, with the Dm ordered in different ways 
(Fig. I, Ref. [^); but (b) only regime I is observed 
in a wide variety of DNSs [2^ l23| . Regime I has 
sufficiently depleted nonlinearity so that a global 
attractor exists, provided the solutions remain in 
this region, as they do in all the DNSs examined 
so far from this point of view. 

The analog of the above theoretical framework 
is developed for the case of 3D MHD turbulence. 
Then the behaviors of D^ are examined - the 3D 
MHD counterparts of Dm in 


21 l23| - 


in a va¬ 
riety of DNSs, which have been carried out inde¬ 
pendently by different groups, to obtain new in¬ 
sights into the depletion of nonlinearity here. It 
is found that 3D MHD turbulence is like its fluid- 
turbulence counterpart inasmuch as all solutions 
remain in only one regime, with depleted nonlinear¬ 
ity, in a large variety of DNSs. The implications 
of our results are also examined for the exponents 
that characterize the power-law inertial range 
dependencies of the energy spectra £^{k) of the 
Elsasser variables on the wave number k. 

The remainder of this paper is organized as fol¬ 
lows. In m the MHD equations are introduced 
and our numerical methods are summarized. 31111 
contains the mathematical analysis of these equa¬ 
tions. 3IVI is devoted to the energy spectra that 
emerge from these calculations. m contains the 
principal conclusions of the paper. The technical 
details of some of our calculations have been rele¬ 
gated to Appendices A, B and C. 


II. MODEL AND NUMERICAL METHODS 


virj. The two forcing functions, f^{x), are de¬ 
fined by 

/±(a^) = A±/,. (5) 

which are absent in decaying MHD turbulence, 
j = V X b is the current density. The mean mag¬ 
netic field Bq in zero in our simulations. The fol¬ 
lowing notation will be used for spatial and tem¬ 
poral averages: 

{■)v=L-^ f-dV, (6) 

Jv 


(•)^ = T-i/' -dt, 

Jo 

with the A^-spatial norm represented by 



( 7 ) 

( 8 ) 


The Taylor-microscale Reynolds number is de¬ 
fined as 


Ra — '^rms^ 




( 9 ) 


with Urms the root-mean-square velocity. Our 
DNSs of the 3D MHD equations use a periodic cu¬ 
bic box and a pseudo-spectral method |l2Hlfil | with 
large-scale initial conditions, and in some cases, 
forcing (Table H] ). All our numerical simulations 
are fully de-aliased. 

For ideal 3D MHD (i.e., i'± = 0, = 0) the 

invariants are the energies 


E± = ^{z±-z±)y=ET±Hc (10) 


A. The equations in Elsasser variables 


together with the magnetic and cross-helicities 


The velocity u and magnetic held b can be com¬ 
bined into the Elsasser variables 

z^ = u±b. (3) 

Then the incompressible 3D MHD equations are 

{dt-^ ■V)z^ = v+V'^z^ +iy^'i7^z^ 

- VV + f^, (4) 

where V ■ z^ = 0, V is the total pressure, i'± = 
i(^ ± 77), and V and 77 are, respectively, the kine¬ 
matic viscosity and the magnetic diffusivity, whose 
ratio yields the magnetic Prandtl number Pm = 


HM = {A-b)y, Hc = {u-b)y, (11) 

where the vector potential A is related to b by 

b = VxA, ( 12 ) 

and the total energy is 

Et = i (m • M -I- 6 • b)y 

= Eu E\,. (13) 

The relative rates of magnetic and cross-helicity 
are also dehned as 

Om = cos{A, b) , ac = cos{u, b) , (14) 
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with |(Tc,m| < 1- These represent the degree to 
which the fields are aligned and they are also mea¬ 
sures, global or point-wise, of the strength of non- 
linearities in MHD. 

By defining the two combinations of the vorticity 
and the current as 




(15) 


it is shown in Appendix that evolve accord¬ 
ing to (with Pm = 1) 


{dt + - lyAu:^ (16) 

3 

= X -|- ^ diz"^ x diZ^ -|- V x . 

i=l 

The two terms on the right-hand side stem from 
the equation for the current; the labels i = 1, 2, 
and 3 refer, respectively, to x, y, and 2 . 

In the ideal case, the constraints that follow 
from conservation laws involve mixed [u, b) cor¬ 
relators 2J, [23 • In the absence of a strong uni¬ 
form magnetic field Bq, magnetic fluctuations at a 
scale comparable to that of the system, Bl, play 
a role equivalent to that of Bq for the small scales, 
provided there is sufficient scale separation, i.e., for 
high-Reynolds-number flows. It has been argued in 
14| that measurable anisotropy develops for scales 
smaller than the Taylor scale based on Bl- There¬ 
fore, the inertial-range energy spectrum can be of 
either Kolmogorov (K41) or Iroshnikov-Kraichnan 
(IK) forms, depending on the cross-correlation. Di¬ 
mensional analysis gives 


C„ = p/4, 


(17) 


if the model of Iroshnikov and Kraichnan (IK) [26 
27| is used and crc = 0; or 


16, 1 


Cp^"'=p/3, 

1^ is used. Appendix [Eld: 


(18) 


iscusses some 


if K41 

of these scaling arguments in a phenomenonogical 


J 


manner. Moreover, 


£^{k) 


fc-3/2 

k-5/3 


(IK) 
(KAl). 


(19) 


Some models and DNS results SSlill 

indicate that the departure from linear scaling, be 
it of the IK or K41 forms, is stronger in 3D MHD 
turbulence than in 3D Navier-Stokes (NS) turbu¬ 
lence, which suggests a depletion of nonlinearity 
by virtue of the tendency of alignment or anti¬ 
alignment of u and b 32 


B. Description of runs 

Table U contains the parameters for the runs an¬ 
alyzed in this paper. All runs have been performed 
in three dimensions by using periodic boundary 
conditions, no imposed external magnetic field and 
a magnetic Prandtl number Pm of unity, except for 
the pm-runs ; no modeling of the small scales is em¬ 
ployed. For the sd-runs (spin-down) the Reynolds 
number is varied. The initial condition for the 
spin-down runs sd is the three-dimensional Orszag- 
Tang vortex [^, with added phase shifts to set 
ac — —0.21 initially. The Aa-Ae runs are high- 
resolution forced runs 35, 3g, with a constant ve¬ 


locity and magnetic forcing for which all the modes 
in the first two Fourier shells are kept constant. 
From the Aa to Ae runs, the resolution increases 
with the Reynolds number. The tg-runs are forced 
in both the velocity and induction equations. In 
these runs, the four-fold symmetries of the Taylor- 
Green vortex extended to MHD are implemented. 
Moreover, the three runs have different resulting 
energy spectra (IK, K41, and k~^), although they 
have the same ideal invariants but with different 
cross-correlations [s^. Finally, the pm runs have 
a fixed viscosity, but variable magnetic diffusivi- 
ties and thus allow for extending the analysis to 
the case of Pm ^ 1 (see |16l|l. 


r 


III. MATHEMATICAL ANALYSIS 


The generalization of the analysis of Refs. [22, 
23| for the 3D NS equations is now described in 
the case of the 3D MHD equations. The relevant 
partial differential equations (PDFs) in Elsasser 


variables are o and (HSl). Two spatially and tem¬ 
porally averaged velocities, C/^, based on are 
defined as 


= L-^ ■ 


ll^>T 


( 20 ) 
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TABLE I: Parameters for our direct numerical simulations, fcmax = N/?> is the maximum resolved wavenumber at 
grid resolution N (the standard 2/3 de-aliasing rule). At and Ra are defined in ([9]). gc and gm are the relative 
rates of cross-helicity and magnetic helicity, respectively. are the parameters extracted from the data for high 
m (subscript max) and low m (subscript min) (See FigO column 3). 


Run 

N 

Ra 

At 

Pm 

n-c 

Gm 

y-l- 

^min 

\ + 

'^min 

"^max 

sdl 

128 

14 

0.27 

1 

-0.27 

-0.22 

1.096 

1.158 

1.101 

1.169 

sd2 

256 

21 

0.20 

1 

-0.27 

-0.23 

1.103 

1.165 

1.116 

1.186 

sd3 

512 

30 

0.15 

1 

-0.27 

-0.24 

1.111 

1.171 

1.129 

1.197 

sd4 

768 

45 

0.11 

1 

-0.26 

-0.24 

1.121 

1.184 

1.141 

1.206 

Aa 

512 

35 

0.098 

1 

0.019 

0.003 

1.049 

1.150 

1.049 

1.156 

Ab 

1024 

54 

0.074 

1 

0.017 

0.004 

1.057 

1.197 

1.060 

1.195 

Ac 

2048 

120 

0.036 

1 

0.011 

Data Not Available 

1.076 

1.167 

1.076 

1.176 

Ad 

2048 

161 

0.027 

1 

0.009 

Data Not Available 

1.074 

1.168 

1.073 

1.157 

Ae 

4096 

341 

0.014 

1 

0.010 

Data Not Available 

1.070 

1.163 

1.072 

1.174 

tgi 

1024 

100 

0.066 

1 

0 

0 

1.121 

1.196 

1.117 

1.197 

tga 

1024 

83 

0.084 

1 

0 

0 

1.161 

1.202 

1.138 

1.202 

tgc 

1024 

110 

0.056 

1 

~ 0.05 

0 

1.084 

1.183 

1.089 

1.175 

pmOl 

512 

240 

0.14 

0.1 

0.122 

0.0047 

1.078 

1.238 

1.078 

1.234 

pm02 

512 

140 

0.10 

1.0 

0.075 

0.0049 

1.070 

1.171 

1.069 

1.160 

pm03 

512 

80 

0.06 

10 

0.226 

0.0077 

1.053 

1.149 

1.052 

1.158 


In turn, the allow us to define two Reynolds 
numbers 


By Holder’s inequality, the are naturally or¬ 
dered such that 


Re± = LU^/u. (21) 

The Reynolds numbers are based on average veloc¬ 
ities, while two Grashof numbers Gr± are based on 
the forcing functions f^{x): 

Gr± =L3/2 ||j± 112/^2^ (22) 


^ ^m+1 ■ ( 25 ) 

If a signal has no intermittency, then the will 
be packed close together, whereas a strongly inter¬ 
mittent signal will cause them to spread out widely. 
The following scaling was first introduced in work 
on the 3D Navier-Stokes equations 2ll-l23l| and will 
be followed here: 


For the class of forcing functions spectrally concen¬ 
trated around a single length-scale {(. = L for the 
purposes of this paper), a relation exists between 
Gr± and Re± for solutions of® derived through 
the method of Doering-Foias [38j (see Appendix ICl) 
where it has been shown that 


where the exponent am is defined as 

2m 


(26) 


(27) 


Gr± < ci?e± (i?e:F-b 1) , Gr± > 1. (23) 

The main variables used in this paper are 
norms of the vorticity field, defined in such a way 
that each has the dimension of a frequency: 

/ r \ 1/2"* 

0±(t)= . (24) 


and where vuo = vL~^ is the box frequency. The 
-scaling comes from symmetry considerations. 
The ordering of the in (|25ll does not necessarily 
hold for the as am is decreasing with respect 
to m. The are the main variables to be used. 
Under the assumption that ® has a solution we 
now look at the evolution of Df : 
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[ \Vuj^\^dV + L-^Wo^ [ \u:^ ■ ■ Vz^)\ dV 

Jv Jv 

3 /. 

‘ V / • {d^z^ X d,z^) I dV + Gr±Df^^'^ . 

.■_i JV 


-3„-3 ' 


+ L w, 


To estimate the first nonlinear term in (I28I) we write (1 < to < oo) 

f f 2m-3 1 2m-3 1 

/ \u}^ ■ {u}"^ ■Vz^)\ dV < / |w=F| 2(^-1) |y_ 2±| dy 


(28) 


< Cl. 


iV 




^\^dV 


2m — 3 
4(m-l) 


dV 


4m{m — 1) 


(29) 


2m — 3 


|a;=F|2m 


4m,(m-l) 




ih 1 2m 


dy 1 


^ l/2m 


Note that the sum of the five exponents in the latter expression is unity. For the last term we have 
invoked the inequality [i^, which requires a Riesz transform in its proof, namely, 


|V2:=‘=||2m < C2,m||‘4;=‘=||2m : 


(30) 


provided 1 < to < oo. Then dMl) becomes 


L ^ 


a; 


im — i 

• V2±)| dV < C3.m [Dt] [Di] 


2cm(m-l) 


2m — 3 

< C3,m[Dt]^^ [D^y 
X [Df] [Df] ^ 


2m-l 

2a^(m-l) 


(31) 


As in the 3D Navier-Stokes equations, this estimate of the nonlinearity is too strong for the dissipation 
terms. However, what was observed in computations of the 3D Navier-Stokes equations is that it displays 
numerically much weaker behavior than the estimate equivalent to m 0- This can be measured by 
numerically tracking Dm in terms of Di, the equivalent of which for 3D MHD isIH^ : 

iA± , 


where, for 2 < to < 9, A* ^ is defined as 


D^<[Diy 

mX^ -b 1 - A=^ 


4± — 

^m,X — 


(32) 


(33) 


4to — 3 

In effect, is a fitti ng p arameter for the maxima in time. An explanation why such a relation should 
hold can be found in [39|. The range of values of A^ have been determined numerically (see Fig. [2] and 
Table m. By inserting (15^ into (1311) it is found that 




X±(2m-l) + m(2m-3) x^ + m(2m-3) 

|a;=^ • (a;=F • Vz^) \ dV < C4,m [Df] . (34) 


Next, the second nonlinear term in (l28l) is considered where (l30l) is used. From this, it is found that the 
estimate for the right-hand side of (IM)) is the same as in (l2^ . apart from the constant C 5 ,m : 


^ 2m-3 

J \u}^ ■ {d^z^ X diZ^)\ dV < C5^m(^j (^J 


4m{m — 1) 


( 35 ) 




=F|2 


2m — 3 

4(m-l) 


O) 


=p \2m 


4m{m — 1) 




± \2m 


l/2m 
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Converting this into the gives the same right-hand side as in (IMl) but with a constant C 2 ,m- Taking 
all these terms together, (E51) becomes 

f X±(2m-l) + m(2m-3) xS+m(2m-3) . 

/ \Vu}^\^dV -|-C6,m 

Jv 

To handle the coupled nature of the ±-variables we define 


^£,±j 4™,™-!, [dT] (36) 


X = Df and Eq = max (37) 

and the two bounded dimensionless energies are defined by \z^\^ dV. By adding the 

±-equations and using the depletion formulas (15^ and (I55)) . a differential inequality is found for X{t) 

+ 2 max (Gr+, Gr_) X^/^ . (38) 

ZhQ 

Note that when = A, as in the Navier-Stokes case, then the exponent of the nonlinear term 

reduces to 1 -I- |A, as it should. 


Without the use of the numerically observed de¬ 
pletion in (1321) , standard methods in analysis leads 
to a term oc X^ in Eg. (1551) (see Ref. for the NS 
case), which does not lead to a control over the so¬ 
lutions at large times. However, provided A^ and 
A::p satisfy 


The e^-term has been left off the figure as it is 
small and can take either sign. When 

2 -h < A=^ < 4 , (regime II) (42) 

there is depletion, but not enough to control solu¬ 
tions; and, finally, when 


1 + < 2, (39) 

4m 

an ‘absorbing ball’ for X exists because Eq is 
bounded above. This ball has finite radius (de¬ 
pending on the upper bound on Eq) into which so¬ 
lutions are drawn if initial conditions are set out¬ 
side the ball, and which cannot escape if initial 
conditions are set inside. Expression (1391) can be 
rewritten as 

A± < 2 + e±, (40) 

where = (A^ — A’*‘)/2m, which is a small num¬ 
ber. Subject to the constraints on A* in (liOl) . the 
ball is such that the 2 ^ are L^-bounded, and thus 
so are u and b. Additionally, the control of X that 
(l38l) affords (an i7i-bound) is also enough to prove 
its compactness. This ball is thus the global at¬ 
tractor which governs the long-time dynamics of 
the PDEs. 

The natural, 3D-MHD analogs of Fig. 1 in 
Ref. [ 2 ^ are the schematic plots of versus 
and iA“ versus in Fig. [U which show three 
regimes. For regular solutions, we must have 

1 < A^ < 2 -I- e* , (regime I). (41) 


> 4 , (regime III) 


(43) 


then > CmDf ■ However, any initial data set in 
this region would be pathological as it would have 
to be prepared as a very large spike in in which 
the L°“-norm is much larger than the L^-norm. In 
the NS-case, it can be shown that solutions are reg¬ 
ular in regime HI, but no more than algebraically 
increasing because of the forcing (39| . Our DNS 


data indicate that regime I (1 < < 2 -|- is 

obtained in 3D MHD for all the solutions we have 
studied. 

In Fig. [^representative results from four of our 
DNSs are given. The first column of Fig. [2] con¬ 
tains log-log (base 10) plots of the energy spectra 
E{k) versus k. Most of these energy spectra show 
power-law forms in the inertial range with an ex¬ 
ponent that is consistent with the K41 value —5/3. 
However, this exponent is consistent with the IK 
value —3/2 for run-tgc; and it is —2 for run-tgi. 
We find that these exponents can depend on the 
values of ac and ctm , which are given in Table |T| 

The second column of Fig. jSjhas plots of A+(t) 
versus t, from which ^ is determined (plots 
for (t) versus t are similar), which follow from 
D+(t). The region where the data lie do not fol¬ 
low exactly the contour boundary curves of Fig. [TJ 







7 



able up to m = 10; note the ordering of the ^ 
is the same for all our runs. 

In the third column of Fig. [21 plots of versus 
m are given, in the range 2 < to < 9, for which 
good-quality numerical data have been obtained. 
From these plots has been found from the min¬ 
imum over TO of A^. In general, 1 < A± < 4; 
however, in all our DNSs, 1 < A± < 2, i.e., our 
solutions lie in regime I. 


FIG. 1: (Color online) Schematic plots of versus 
showing the three regimes in 3D MHD (see text). 
The additive term in (SB has been omitted because 
it is small and can take either sign. The values of A* 
in Table [I] lie only just above the lower bound A^ = 1. 
Solutions are regular in regime I but not in regime II. 
To start in regime III requires unphysical initial data. 


IV. SPECTRA 

A. How to estimate the spectrum for the 3D 
MHD-Elsasser system 


The A^ have been determined as in Ref. for 
the 3D-NS equations. are defined to be those 
values that have been computed from Eq. (1551) for 
A ■ ^ check has shown that our data are reli- 


The method of Doering and Gibbon is now 
followed which explains how to estimate average 
length scales and a corresponding spectrum based 
on ideas in [I|. It is necessary to define a set of 
time-averaged inverse length scales |51| 


/.2 ±2x /L^Vu;mdV\ / 

■ \ ll-^lli /t-\ ^ , 


(44) 


where the labeling of the subscripts is based on the number of derivatives on z^. Dividing (1561) by 
and time averaging, we find 521 



^ ^ x^,+ml2r,^-3) \ 


< ce,m ( [Dt] 

[Df] \ 

/ T 



xS+m(2m-3) 


< Ce,ra{Dt)T 


(45) 


Secondly, it is necessary to estimate by us¬ 

ing the energy inequality version of (|3]) 

^^||2:=^ll2 < -I^ll‘^'^ll2 + I|2:=^I|2||/'^||2 ■ (46) 

By time averaging, converting into a dimensionless 
form, and using (1551) . it is found that 

{D^)j,<Gr±Re±<cRel{Re^ + l) . (47) 

Moreover, by introducing the definitions, the first 
of which is the Elsasser analog of the Taylor micro¬ 


scale. 


^±2 _ 
^ 1,0 ~ 


±4 l|Va;±||i 

IU±I|2 ’ ^^2,0 |U±||2 ’ 


and adapting ideas in [^, (1551) gives 

Then it is easily shown that 


(48) 


(49) 


( 50 ) 























FIG. 2: (Color online) Representative results from our DNSs: total energy spectra (first column), temporal 
evolution of (second column) and values of A* (third column). The rows correspond to different runs: 
a) decaying 3D MHD turbulence; b) forced, statistically steady 3D MHD turbulence; c) forced, statistically 
steady 3D MHD turbulence with imposed Taylor-Green symmetries and d) forced, statistically steady 3D MHD 
turbulence with Pm =0.1. For parameters see Table HI 


and so from (l44l) . (|4^ and (ITZl) . in which only the 
dominant term has been kept, it is found that 

< C6.mi?e±™ (i?e:p -b 1)^™ , (51) 

where, with xtk = {'ni — 1)A^ -I-1, 


— 


Pm 


Thus, 


A=^-bl /A=F-A± 

4 \ 8m 

(2m - 1)A=F -b A± 
16m 

can be written as 

< cRe^ “ Re 




(52) 

(53) 

(54) 


where the factor of unity has been ignored in the 
large i?e^-limit and the limit of large m has been 
taken. Then the problem is whether it is possible 
to consider Re+ and Pe_ as independent variables 
or not. The simplest way is to note that 

< ||v||2 + 2||v||2||b||2 + ||b||l 

< 2 (11^111 +||6||^) =4Ptot. (55) 

Defining a global Reynolds number as 


Re = L^j2Etotlv , 


(56) 
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- m=1 

- m=2 

- m=3 

- m=4 

- m=5 

- m=6 

- m=7 

- m=8 

- m=9 

- m=1 

- m=2 

- m=3 

- m=4 

- m=5 

- m=6 

- m=7 

- m=8 

- m=9 


FIG. 3: (Color online) Plots versus Ra of „„ > 

for the runs Aa-Ae (pentagrams) and its decaying- 
MHD analog, for the runs sdl-sd4 (circles), where we 
use the value of at the time at which the energy 
dissipation rate e reaches its hrst maximum [l^. See 
Table U for additional information about these runs. 


gives 

/ 4- \ >=*= + ! I AT 

<cRe 4 + 8 . (57) 

Now the implications of these results are examined 

for energy spectra derived in (I57I) . By assuming 
isotropy and the power-law Ansatze 

= (58) 

I 0, k > , 

the identification (see Appendix B and [l|) 

~ (^, (59) 

allows us to find an inequality relation between 
and A±. 


> 3 - 


4 ^5 

2A± -b AT - 3’ 


(60) 


which excludes the IK exponent 3/2, at least in the 
absence of intermittency. However, full isotropy 
has been assumed at all scales and the limit Re —>■ 
oo when comparing and (I551) . Furthermore, 
(ICT)) implicitly assumes E+ = E-. In general, this 
relation is modified leading to a set of inequali¬ 
ties for the q^ that do not exclude IK (the last 
two equations in Appendix B). In fact, IK is ex¬ 
cluded only if correlation is neglected. This is con¬ 
sistent with the fact that run-tgc (see Fig. 2 and 
Table I) produces an IK scaling because it has a 
non-zero cross-correlation. It is known that in the 
presence of cross-correlations between the velocity 
and magnetic field {Elc ^ 0), different indices arise 
for the energy spectra of the two Elsasser fields, 


(see, e.g., 


4l| and references therein), a result 


that persists in the case of weak MHD turbulence, 
as shown through wave-turbulence developments 


42| |. The mathematical analysis as well as the nu¬ 


merical simulations presented in this paper, all put 
on firm ground that indeed Hq plays a crucial role 
in determining the distribution of energy among 
scales. 

In Fig.[3]plots of versus Ra have been 

displayed, where the angular brackets now indi¬ 
cate the average over the statistically stationary 
turbulent state; we present data (pentagrams in 
Fig. from the runs Aa-Ae (Table |T|. Circles in¬ 
dicate decaying-MHD data points (runs sdl-sd4); 
here we use the value of at the time at which 
the energy-dissipation rate e reaches its first max¬ 
imum |l6l| . The order-m moments of gradients, 
or gradmoments, of hydrodynamic fields have been 
used to investigate Nelkin scaling ^-4^, i.e., the 
power-law dependence of the gradmoments on the 
Reynolds number Re in the case of fluid turbu¬ 
lence; the Nelkin-scaling exponents can be re¬ 
lated to the structure-function exponents Crr 


43- 


45l| . Our vorticity moments are upper bounds for 


gradmoments of the Elsasser variables (see (1501) ). 
If these bounds are saturated, then the exponents, 
which can be extracted from the plots of Fig. [31 
should be related to the Nelkin exponents for 3D 
MHD turbulence. A more detailed exposition of 
such scaling in 3D MHD has been deferred to an¬ 
other study. 

In liquid metals, as well as in the solar photo¬ 
sphere, Pm is very small. It can also be very large 
as, e.g., in the interstellar medium. Our mathe¬ 
matical analysis is not valid if Pm ^ I because v- 
can become negative for Pm < 1. However, our 
DNS results in Fig.|3]show that plots for Pm = 0.1 
(bottom row) are similar to their counterparts for 
Pm = 1 (top three rows). Furthermore, at least for 
a hxed value of ly, Table U shows that the values of 
are comparable to their Pm = 1 counterparts; 
these values of A^ decrease marginally as Pm is 
increased. 


V. CONCLUSIONS 


Our work, which builds upon the studies of 
Refs. 22, [ 2 ^ for fluid turbulence, provides 
insights into the depletion of nonlinearity in 3D 
MHD turbulence and its intermittency. In par¬ 
ticular, we have introduced the scaled moments 
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, and then obtained inequalities containing 
and Df ; these inequalities specify three possible 
regimes. In essence, it has been found that 3D 
MHD turbulence is similar to its fluid-turbulence 
counterpart insofar as all solutions that have been 
investigated have remained in only one regime 
(regime I), which displays depleted nonlinearity 
(Fig. [1]). Moreover, under the assumption of 
isotropy our results lead to the inequality (IM)) for 
the spectral exponents q^. In fact, the inequal¬ 
ity (1301) can relate and the order-m moments 
of gradients of the magnetohydrodynamic fields; 
such moments can then be used, along with a suit¬ 
able generalization of Nelkin scaling 43l-l45l| for 3D 
MHD turbulence, to relate slopes of plots like those 
in Fig. [3] to the multiscaling exponents of Elsasser- 
field structure functions. We conclude that 3D 
MHD appears to have more nonlinear depletion 
than fluid turbulence, because the values of are 
lower than those for their fluid-turbulence coun¬ 
terparts ; this can be attributed to Alfven waves 
weakening the nonlinear eddies. 


Acknowledgments 

We thank H. Homann and R. Grauer for pro¬ 
viding the data of runs Ad and Ae. Data from 
runs Ad and Ae of H. Homann and R. Grauer 
were produced on the IBM BlueGene/P computer 

I 


JUGENE at FZ Jiilich made available through 
the “XXL project of HB028”. We also acknowl¬ 
edge U. Frisch and D. Vincenzi for useful discus¬ 
sions. The computations for runs sdl-sd4 were 
performed on Janus (UC-Boulder). Runs Aa-Ac 
and tg were performed at Mesocentre SIGAMM 
hosted at the Observatoire de la Gote d’Azur and 
GIGADA hosted by the University of Nice-Sophia. 
Gomputer time was also provided by GENCI on 
the IDRIS/GINES/TGGG clusters. JG and AP 
thank Eederation Doeblin for support. RP thanks 
the Department of Science and Technology (India) 
for support and SERG (HSc) for computational re¬ 
sources. AG is grateful for support through a grant 
from the European Research Council (ERG) under 
the European Community’s Seventh Framework 
Programme (FP7/2007-2013)/ERC Grant Agree¬ 
ment No. 297004. GS acknowledges support 
from the ERG Advanced Grant “NewTURB”, No. 
339032. GK thanks the Indo-French Centre for 
Applied Mathematics (IFCAM) for supporting a 
visit during which parts of this paper were writ¬ 
ten. JS is supported by the National Science Foun¬ 
dation (NSF) Graduate Research Fellowship Pro¬ 
gram (GRFP) under Grant No. DGE 1144083. 
AP and JDG acknowledge, with thanks, IPAM 
UCLA where this collaboration began in the Au¬ 
tumn of 2014 on the program “Mathematics of 
Turbulence”. 


Appendix A: Vorticity and current equations 


The following four identities have been used for vectors D and Q (not necessarily divergence-free): 

V X [V X D] = V[V • D] - 

V[JD -Q] = Dx[VxQ] + Qx[VxD] 

-t D VQ -b Q VD 
V[£)xQ] = QS/xD-DVxQ 

\7 x[DxQ] = Q VD - D SJQ . (Al) 

The equation for the vorticity is derived straightforwardly: 

{dt + u- V)u; = us ■ Vu -(- b ■ Vj — j ■ Vb. (A2) 


By using D = u x b in the above identities (with V • D 7 ^ 0), we obtain the equation for the current 

dtj = V [V • [it X 6] - V^[it X b]] , (A3) 

which, upon expansion, leads to : 

V[V • [it X b]] = uj ■ Vb -f b • Vu; — b x V^it + uj x j — u ■ Vj — j ■ Vit -|- it x V^b — j x u) (A4) 
— V^[it X b] = —V^uxb — uxV^b — 211idiUxdib. (A5) 
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So the equation for the current (note the cancellations in the terms) is: 

{dt + u ■ V)j — (jJ ■ V6 — b ■ Vo; + j ■ Vit — 2a; x j = —2Y,idiU x dib , 


(A6) 


an expression already written in For the curl of the Elsasser field a;+ = uj + j, (with ± symmetry 
for dtu}~), (IA6I) and (IA2I) reduce to (ITCl) 


{dt 


V)a;+ = a; • 


u) XU}' 


T^idiZ^ X diZ 


(A7) 


Note that the geometry term 2ujxj in the equation 
for the current density does not appear in the vor- 
ticity equation; also, it is weak for almost-aligned 
current and vorticity (or o;^) The equa¬ 

tions (TATI) for the temporal evolution of u)^ fol¬ 
low immediately from the above. Note also that 
a;~ X a;+ = 2a; x j does not affect the point-wise 
production of o;^, whereas the second term can 
create current density; here the labels i = 1, 2, 
and 3 refer respectively to the x, y, and z deriva¬ 
tives. Finally note that, for a flow evolving towards 
strong local correlations between the velocity and 
magnetic field (z+ = 0 or 2 :“ =0), this extra term 
is weak. 


Appendix B: Phenomenological argument for 
fluids and MHD 

In the fluid case, the total energy and dissipation 
can be written in terms of the energy spectrum 
with spectral index q as 

= f Ak-^ , e = z/ / ° , (Bl) 

J kQ J ko 

with the dimension of A as [A] = [e“][A^]. One 
finds straightforwardly a = 2/3, 6 = [5 — 
3 ( 7 ]/ 3 , so & = 0 for (7 = 5/3, as expected. 
This leads to a cut-off wavenumber fc^/fco = 
fcg or in terms of the 

Reynolds number, 

Re = UL/v = Lkc = Re^, 

X = [3-(7]-\ (B2) 

with ko = 27r/L. 

In MHD, one can follow the weak-turbulence IK 
prescription (remaining in the isotropic framework 
for simplicity). Then, A = where Bq is a 

large-scale strong (quasi)-uniform magnetic field; 
so c = 1/2, d = [3 — 2q]/2. This leads to 

kjko = fc-3/[2(3-<?)] , 


or, in terms of the Reynolds number. 

Re = UL/v = \eBQL^\^l^v -^, (B4) 


with 


Lkc = Re^, a: = [3 — g] 

r = ^ « I , (B5) 

Bq 

which is a hypothesis that is compatible with the 
wave-turbulence assumption. Thus, with the in¬ 
troduction of the factor r, the scale dependence of 
the cut-off wavenumber with Reynolds number is 
the same for fluids and MHD. 

This phenomenological argument can be repro¬ 
duced in the more general case when the veloc¬ 
ity and magnetic fields are correlated, i.e., with 
A+ 7 ^ E-. This results in a condition between the 
indices q± of the E± spectra, namely g+ -I- g_ = 3 
within the IK framework with, as before for the un¬ 
correlated flows, q+ = q- =3/2 (see [i^ for an in¬ 
troduction). Two-point closure computations and 
two-dimensional numerical simulations find q+ 7 ^ 
g_ at high correlations, but the three-dimensional 
case remains open. As in the preceding case of un¬ 
corrected MHD, we make the assumption that the 
± integral scales are both comparable to the box 
size L. 

After some algebra along the same lines as be¬ 
fore, one finds that the dissipative wave numbers 
for the E±{k) spectra are equal both to k^ = k_ = 
k, = e/[z;2Ho]i/3, as in the zero-correlation case, 
or: 

Lk, = [z+/Bo]^/^[zo/Bo]^/^Re]/^Re]/^ , (B 6 ) 

with Re± = z^L/v, and we have assumed that the 
magnetic Prandtl number is equal to unity so that 
V = rj = = v_. Writing 

E±{k) = A±(e,Bo, T)fc“''±, 
A±(e,Ho, A) = 


(B7) 
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it is readily found that, under the assumption that 
6+ + 5_ =-3(a+ + a_)+4=l (B8) 

(so that E+{k)E-{k) ^ [eBo]k~^), and that 

C-^. + c_ = —3 + (3 — <7+) (3 — Q—), (B9) 

then 

Lke ^ ) + ^ (BIO) 

or, in terms of the Reynolds numbers Re±, 

Lk, ~ [[z+/Bo][z-/Bo]]^^^^ETi^ 


X [i?e+i?e_]. 


(Bll) 


Finally, this phenomenological relation can be used 
as in the non-correlated case to establish the upper 
bounds of spectral indices. In this case one hnds 
that 


[1-1(9+-!)] <^i^rA+ + l A_ 


(3 - q+) + (3 - q-) 
and similarly 


< min 


A— + 1 A-i 


(3-9+) + (3-9-) - I 4 

(B13) 

with 1 < 9 + < 3 and 1 + A±/2 < 2 using equa¬ 
tion (|40)) . These results are to be contrasted with 
the uncorrelated case obtained in the Kolmogorov 
framework. It is possible to show that in the cor¬ 
related case the IK spectrum 9 + = 3/2 cannot be 
excluded. 


Appendix C: The Doering-Foias Gr± - Re^ 
relation for MHD 


Following Doering and Foias 


the forcing 


function f^{x) is split into its magnitude E^ and 


its “shape” (p^ such that 

f^{x) = F^(p^{£-^x), 


(Cl) 


where i is the longest length scale in the force and 
is taken to be .^ = L for convenience in the rest of 
the paper. On the unit torus in d-dimensions, 
</> is a mean-zero, divergence-free vector field with 
the chosen normalization property 


lid 


[v^V^r d‘^y = ^- 


L^-norms of on are 


y^f^\\l = c 


i 0—2Nj^d^-^2 




(C2) 


(C3) 


where the coefficients C^, which refer to the shape 

(C4) 


of the force but not its magnitude, are 


(BI2) Various bounds exist such as (among others) 


||VA-^/±||oo = D±K£ 


The energy dissipation rate e is 


2M-1 


(C5) 


e^ = (^iyL-^J \Vz^\^dv'^ =izL-^{H^) . 

(C 6 ) 

In terms of the Grashof number in be¬ 
comes (f = L) 


Gr± = F^l^/v^ . 


(C7) 


Following the procedure in [3^ (pg 296 equation 
(2.9)), we multiply (|4|) by and inte¬ 

grate to obtain 


^ [ z^-[{-A-^)f^]dV = V j Az^ ■[{-A-^)f^]- 

Jid Jid 

+ j f^-[{-A-^)f^]dV. 


z^ ■ Vz^ ■ [(-A-^)/±] dV 


(C 8 ) 
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Now if we integrate all the terms by parts, and take the time average, we get 


|V-Mj±|2 . < y/L- 




± . (_A-*^+i)/±| dV 


+ \z^ ■[V[{-A-^)]f]-z^\ dv'^ . 


Thus, after a Schwarz inequality, (lea turns into 

By using (IC7I) . in the limit Gr± oo, (IC10|) becomes 

Gr± < c (i?e± + i?e±i?e;p) . 


(C9) 


(CIO) 


(Cll) 
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